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Pair-distribution functions of correlated composite fermions 
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Pair-distribution functions g(r) of Laughlin quasielectrons (composite fermions in their second 
Landau level) are calculated in the fractional quantum Hall states at electron filling factors v e — 4/11 
and 3/8. A shoulder in g(r) is found, supporting the idea of cluster formation. The intra- and inter- 
cluster contributions to g(r) are identified, largely independent of v £ . The average cluster sizes are 
estimated; pairs and triplets of quasielectrons are suggested at v & = 4/11 and 3/8, respectively. 
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I. INTRODUCTION 

Pan et alA have recently observed the fractional 
quantum Hall (FQH) effect in a spin-polarized two- 
dimensional electron gas (2DEG) at the z/ e = |, and 
y| fillings of the lowest Landau level (LL). In the com- 
posite fermion (CF) model, 4,5 these values correspond 
to the fractional fillings v = \ >1 |, and | of the sec- 
ond CF LL, respectively. In Haldane's hierarchy picture 6 
of these states, Laughlin quasielectrons (QE's) fill (the 
same) fraction v of their LL. The most striking conclu- 
sion from Pan's discovery is that the CF's (or QE's) can 
also form incompressible states when partially filling a 
LL. This could not be predicted by a simple analogy 
with known fractional electron liquids (Laughlin, 3 Jain, 4 
or Moore-Read^ states), because of a different form of 
QE-QE interaction, 8 910 therefore yielding qualitatively 
different QE-QE correlations. 

Although several numerical studies of interacting QE's 
have been reportedi ^* 1 ^* 1 ^ and ideas such as CF flavor- 
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mixing,"" QE pairing^* 1 * or stripes^ were invoked, the 
correlations responsible for the FQHE at v e = j- and | 
are not yet understood. It has not even been settled if 
these FQH states are isotropic, and the energies of liquid 
and solid phases were compared recentljiS . (although the 
Laughlin form was arbitrarily assumed for the liquid). 

Sometimes overlooked is a general connectioniSiiS be- 
tween the form of Haldane pseudopotential^ occurrence 
of Laughlin correlations, and the validity of CF trans- 
formation. Actually, the form of QE-QE interaction is 
known from independent calculationst&Sii " and Laughlin 
correlations among the QE's have been ruled out using 
both a general pseudopotential argument 9 and a direct 
analysis of many-QE wavefunctions.— In this paper we 
refer to the following well-established facts: 

(i) The QE-QE Haldane pseudopotentialSi is known 
from exact diagonalization of the Coulomb interaction 
among electrons in the lowest LL&2ii£ Since there are 
no unchecked assumptions in such a calculation, it must 
be regarded a "numerical experiment." Neither finite- 
size errors, lowest-LL restriction, finite 2DEG width, nor 
other details of realistic experimental systems affect the 
dominant feature of the pseudopotential which is the lack 
of strong QE-QE repulsion at short range. 
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(ii) The QE's do no1*iH have Laughlin correlations at 
v = | corresponding to u e = •jp The Moore-Read half- 
filled state is no&2£ an adequate description of QE-QE 
correlations at v = | corresponding to v e = |. 

(iii) A sequence of nondegenerate finite-size QE ground 
states with a gap, extrapolating to v = | has been 
foundiS on a sphere. Although spherical geometry is not 
adequate for studying crystal or other broken-symmetry 
phases, the identified states appear incompressible and 
have the lowest energy of all QE liquids (considerably 
below the Laughlin state). 

To address the problem of correlations at v e 
and we calculate pair-distribution functions g(r) in 
the incompressible liquid ground states of up to N = 14 
QE's. Their comparison with the (known) curves of the 
Laughlin and Moore-Read states implies a different na- 
ture of the QE correlations in these novel FQH states. 
It shows that their incompressibility cannot be explained 
by a simple analogy between the QE and electron liq- 
uids, and suggests that different wavefunctions need be 
proposed for the correlated CF's. Unfortunately, the cal- 
culated g(r) are of little help in a precise definition of 
these wavefunctions, even though some qualitative state- 
ments can be made about the QE correlations. 

From our finite-size results we identify and analyze the 
size-independent features in g(r): the ~ r 2 behavior at 
short range and a shoulder at a medium range, and argue 
that they are consistent with the ideai^ of QE cluster 
formation. Short- and long-range contributions to g{r) 
are found, describing correlations between the QE's from 
the same or different clusters. Both intra- and inter- 
cluster QE-QE correlations depend rather weakly on v. 
The average size of the clusters is estimated; it seems 
that the QE's form pairs at v — | and triplets at v = ~. 
A similar analysis of g(r) carried out for the Moore-Read 
state reveals a qualitatively different behavior. 



II. MODEL 



A. Haldane sphere 



The numerical calculations have been carried out in 
Haldane's spherical geometry^ convenient for the exact 
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study of short-range correlations. In this model, the low- 
est LL for particles of charge q is a degenerate shell of 
angular momentum I = Q. Here 2Q is the strength of 
Dirac monopole in the center of the sphere defined in the 
units of elementary flux 4>q = hc/q as 2Q(j>o — AttR 2 B, 
the total flux of the magnetic field B through the surface 
of radius R. Using the usual definition of the magnetic 
length, A = ^Jhc/qB, this can be written as IX 2 = R 2 . 
In the following, A denotes the QE magnetic length cor- 
responding to the fractional charge q = — e/3. 

The relative (TZ) and total (L) pair angular momenta 
are related via L = 21 — TZ. For fermions, TZ is an odd 
integer, and it increases with increasing average pair sep- 
aration U (r 2 ). The interaction (within the lowest LL) is 
entirely determined by Haldane pseudopotential defined 
as the pair interaction energy V as a function of TZ. 
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FIG. 1: Excitation energy spectra (energy E as a function 
of total angular momentum L; Eo is the ground state energy) 
of N interacting QE's on a sphere, at the values of CF LL 
degeneracy V = 21 + 1 corresponding to the incompressible 
ground states at the QE filling factors v = | (a) and ^ (b). 



B. Exact diagonalization 

Recently, we have identified^ the series of finite-size 
spin-polarized states that in the thermodynamic limit 
describe the FQHE at v e = jj and |. To do so, we 
have carried out extensive exact-diagonalization calcu- 
lations for interacting QE's (particles in the second CF 
LL). On Haldane sphere, N fermionic QE's were con- 
fined in a standard way to an angular momentum shell 
of degeneracy T = 21 + 1, corresponding to the QE filling 
factor v ~ N/Y, and the Haldane QE-QE pseudopoten- 
tial V(1Z) was taken from earlier calculationsA2iiS 

Regardless of the electron layer width w, magnetic field 
B, or other experimental parameters, the dominant fea- 
ture of V(TZ) is strong repulsion at TZ — 3. This feature 
alone determines the wavefunctions at | < v < | (with 
the QE-QE correlations consisting of maximum possible 
avoidance of Haldane pair amplitude Q at TZ = 3), which 
are hence virtually insensitive to the (sample-dependent) 
details of V(TZY This justifies model calculations using 
V(TZ) of Refs. l8Tl9lllCl Actually, a model pseudopotential 
as simple as V = 5-r,3 is sufficient to reproduce correct 
correlations and incompressibility at v e = -A or |. 



III. NUMERICAL RESULTS 

A. Energy spectra 

The numerical results carried out for TV < 14 (two sam- 
ple spectra are displayed in Fig. ^) showed^ a sequence 
of nondegenerate (i.e., at the total angular momentum 
L = 0) ground states at 21 = N / v — 7 with v = i and 
7 = 7. The significant and well-behaved (as a function of 
N) excitation gap along this sequence strongly suggests 
that it represents the infinite v e = jj FQH state observed 
in experiment i The value 7 7^ 3 precludes Laughlin cor- 
relations among QE's in this state (earlier ruled out in- 
directly, based on the form of QE-QE pseudopotential 9 ), 



i.e., an idea that the v e = jj state is simply a Haldane hi- 
erarchy state of Laughlin-correlated CF's. While the ex- 
act correlations in this (known only numerically for a few 
consecutive N) ground state have not yet been defined, 
their vanishing degeneracy (L = 0) implies that they de- 
scribe a QE liquid, rather than a broken-symmetry state 
(such as liquid-crystal nematic states proposed^ in the 
context of FQHE at different values of v) . 

Another sequence was anticipated at 21 = 2N — 7 to 
represent the infinite v e — | FQH state. However, the 
only ground state with a significant gap and remaining 
outside of the v = ^ sequence (or its particle-hole sym- 
metric v — I sequence at 21 — |7V + 2) occursi^ for 
N = 14 and 21 = 25 (and it also has L — 0). These 
values of (N, 21) happen to belong to a 21 = 2N — 3 series 
representing the Moore-Read (pfaffian) paired state, but 
the overlap between the two turns out nearly zeroiiSi 2 ^ 
Moreover, the ground states for the two neighboring even 
(as appropriate for a hypothetically paired state) values 
of N = 12 and 16 (and 21 = 21 and 29) have L > 
and no gap, the value of 21 = 17 for N — 10 coincides 
with the v — | sequence (so that only for N > 8 can 
the filling factor v be meaningfully assigned), and we are 
unable to compute the spectra for N > 18. Nevertheless, 
despite little evidence available from numerical diagonal- 
ization, the ground state for N = 14 and 21 = 25 (and 
its particle-hole counterpart at N — 12 and the same 
21 = 25) may possibly represent the v e = | FQH state 
(i.e., have similar correlations causing incompressibility). 

B. Pair-distribution functions 

The QE-QE pair-distribution functions g(r) have been 
calculated for the incompressible many-QE ground states 
as expectation values of the appropriate pair interaction, 

g(r) = (2/N) 2 (5(R9 - r)) . (1) 

Here, 9 is the relative angle on a sphere, so that r mea- 
sures interparticle distance along the surface (rather than 



3 




2 46 80 2 46 8 

r/X r/X 



FIG. 2: QE-QE pair-distribution functions g(r) of the in- 
compressible ground states at different QE filling factors v. 
(a) curves for v = ~ and different QE numbers N; (b) curves 
for QE's at different v (thick lines) compared to some known 
incompressible states of electrons. 

chord distance). More accurately, r is the distance be- 
tween the centers of extended QE's (note that in the cal- 
culation of many-QE wavefunctions, the system of QE's 
is mapped onto the lowest LL of point charges interact- 
ing through an effective pseudopotential). The prefactor 
in Eq. JQ) ensures proper normalization, <?(oo) — * 1. De- 
noting infinitesimal area by dS = 2irR 2 g?(cos 8) or (in 
magnetic units) by ds = dS/2ir\ 2 , we get an equivalent 
normalization condition, 

J[l-g(r)]d8=^^ V - 1 (2) 

in large systems. Since ds — I d(cos 8), a "local filling fac- 
tor" can also be defined as v(r) — dN/ds — (N/21) g(r), 
and it satisfies u(oo) — v and J v(r)ds = N — 1. 

The results for the v = | sequence at 21 — 3N — 7 
are shown in Fig. |2Ia). Similarity of all four curves 
is evident, indicating size- independent form of correla- 
tions (hence, describing an infinite system), with a well- 
developed shoulder around r ~ 2.5A. Similar shoulders 
occur in g(r) of all incompressible ground states at v = | 
or i (the v = | sequence at 21 = |iV + 2 is obtained 
from 21 = 3N — 7 by replacing N with T — N, while at 
v = i there are two particle-hole conjugate sequences 

at 21 = 2N - 3 and 27V + 1, denoted by v = | ± ). The 

four curves representative of v = |, |, and ^ are shown 
in Fig. Etb). They are all clearly different from those 
marked with thin lines and describing correlations known 
for other incompressible FQH states (full LL, Laughlin 
v = | state, or Moore-Read half-filled state). This is a 
direct indication of a different nature of QE-QE correla- 
tions responsible for the FQHE at v e = jj and |. 

Let us stress that although the QE-QE interactions are 
not known with great accuracy, the correlation functions 
in Fig. |2 are rather insensitive to the details of V(TZ), as 
long as the dominant repulsion occurs at 1Z — 3 (which 
seems to be universally true in the systems studied ex- 
perimentally). This insensitivity is reminiscent of the 
Laughlin wavefunction, which also very accurately de- 




FIG. 3: Gaussian deconvolution of the QE-QE pair distri- 
bution functions g(r): dots - data of Fig. Etb); lines - fits. 



TABLE I: Gaussian deconvolution parameters for QE-QE 
pair-distribution functions shown in Figs.^b) and|3| 
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scribes the actual v — i ground state for a wide class of 
electron-electron pseudopotentials. However, while the 
avoidance of TZ = 1 by the electrons in the lowest LL 
can be elegantly described by flux attachment in the CF 
picture, no similar model has been proposed yet for the 
avoidance of TZ = 3 by the QE's. Therefore, knowing the 
g(r) curves of QE's and understanding their correlations, 
we still cannot write their wavefunctions. 



C. Gaussian deconvolution 

The curves of Fig-Etb) can be accurately deconvoluted 
using gaussians, G(r/\) = Acxp[— (r/X — S) 2 /2a 2 }. This 
is shown in Fig. [3] where the symbols mark the exact 
data of Fig. Hfb) and the lines give the (nearly perfect) 
fits using three gaussians, g = 1 — Go — G\ — G2 (sufficient 
for r < 6A). The fitted values of [ylj,5j,aj] for all four 
curves are listed in Tab.||| Note that Aq — 1, So — 0, and 
61 = 3 for all curves (the latter value being least obvious, 
but probably resulting from the avoidance of the same 
IZ3 = 3 by the QE's at all values of v). The values 
of the G2 parameters are not very meaningful when the 
next term in the approximation (G3) is neglected. The 
clearest difference between the four curves is in A\. 



D. Short/long-range deconvolution 

It appears more physically meaningful to decompose 
g(r) into go = I — cxp(— r 2 /2A 2 ) describing a full lowest 
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FIG. 4: (a) Ratio of QE-QE pair-distribution functions g(r) 
to go(r) of a full lowest LL for different incompressible QE 
ground states; (b) the "remainder" ffdiffO") defined by Eq. 

LL24 and a (properly normalized) "remainder" gas, 

9{r) = ag (r) + (1 - a)g diS (r). (3) 

For each g{r), parameter a is calculated as the limit of 
g/go at r — ► 0. It is clear from Fig. Ufa) that <?(r) is 
accurately approximated by a g (r) within a finite area 
or a radius ~ A for all four ground states of Fig. Hfb). 
The numerical values of a are 0.772, 0.804, 0.856, and 
0.899 for v = |, | , 5 + , and |, respectively. Evidently, 

a is size-dependent (e.g., the pair of values for v = ^ 
must converge to the same thermodynamic limit). 

The four curves gdis(r) calculated from Eq. J3J are 
plotted in Fig. Efb). Symbols and lines mark the exact 
data and the three-gaussian fits of Tab. |]J respectively. 
We note that: (i) For the pairs of particle-hole conjugate 
states (N = 12, 18 at 21 = 29 and N = 12, 14 at 21 = 25), 
the c/diff (r) curves are identical, (ii) The curves obtained 
for v — i and | are very similar (and possibly identical 
in large systems); they all vanish at short range and have 
a minimum at r ~ 3A and a maximum at r w 5.5A. 



IV. DISCUSSION 

A. QE clustering 

Some information about the form of QE-QE correla- 
tions can be easily deduced from the form of interaction 
pseudopotential V(TZ), which is simply the interaction 
hamiltonian defined only for those pair states allowed in 
the lowest LL. In low-energy many-body states the par- 
ticles generally tend to avoid pair eigenstates with high 
interaction energy, which means minimization of the cor- 
responding Haldane pair amplitude Q. If repulsion V de- 
creases sufficiently quicklj*i& as a function of TZ (the exact 
criterion being 19 that V decreases sublinearly as a func- 
tion of y/ (r 2 )), the smallest value of TZ — 1 is avoided. 
This Laughlin type of correlation is elegantly described 
by attachment of 2p — 2 fluxes to each particle in the 
CF transformation. In a Laughlin-correlated state, each 



particle avoids being close to any other particle (as much 
as possible at a given finite v). 

When short-range repulsion weakens (V at TZ = 1 de- 
creases compared to V at TZ > 3), Laughlin correlations 
disappear and can be replaced by pairing or formation 
larger clusters. Pairsi^i^ or clustering^ were suggested 
by several authors for the QE's. This idea was justified by 
an observation that QE-QE pseudopotential nearly van- 
ishes at TZ — 1 and is strongly repulsive at TZ = 3, causing 
an increase of Q(l) and a simultaneous is decrease of Q(3) 
compared to the Laughlin-correlated stated 

The assumption that QE's form clusters naturally ex- 
plains a shoulder in g(r), and allows one to interpret 
go and <?diff as the intra- and inter-cluster QE-QE cor- 
relations, i.e. the short- and long-range contributions 
to g, corresponding to the QE pairs belonging to the 
same or different clusters, respectively. The vanishing of 
Staff ( r ) at short-range reflects isolation of QE's belonging 
to different clusters. The reason why g^is is not positive 
definite is that intra-cluster correlations are accurately 
described by go only within a certain radius. In other 
words, the actual inter-cluster contribution to g is not ex- 
actly given by <7diff defined by Eq. 10 . Nevertheless, the 
following two conclusions remain valid: (i) the intra- and 
inter-cluster QE-QE correlations are similar at v = -| , | , 
and |, with the respective correlation-hole radii go ~ A 
and Q\ ~ 4A; and (ii) the cluster size K depends on v. 

Similar form of g{r) was found 23 for broken-symmetry 
Laughlin states, in which the shoulder results from an- 
gular averaging of an anisotropic function g(r, tfi) ~ r 2 or 
r 6 , depending on <f>. However, the present case of QE's 
is different, because g(r) is isotropic (wavefunctions have 
L = 0) and the shoulders result from radial averaging of 
inter- and intra-cluster correlations, (beginning as ~ r 2 
and a higher power of r at short range, respectively). 

B. Average cluster size 

In a clustered state, the (average) cluster size K is 
connected to a, and the form of <?diff depends on correla- 
tions between the clusters. The values of K at v — i or i 
can be estimated by comparison of the actual parameters 
a with those predicted for the hypothetical states of N 
particles arranged into N/K independent i^-clusters. By 
independence of the clusters we mean that inter-cluster 
correlations do not affect the local filling factor v{r) at 
short range. For a single cluster, which on a sphere is the 
if-particle state with the maximum total angular mo- 
mentum L = Kl — \K{K — 1), the vk{t) depends on the 
surface curvature and thus (through R/\ = y/l) on 21. 

We have calculated the prefactors (3k of the short- 
range approximation z^R-(r) « /3x5o( r ) f° r different val- 
ues of K and 21 and listed some in Tab. |n] (note that 
v-i(r) is known exactly). These coefficients are to be 
compared with (3 — (N/21) a of the incompressible N- 
QE states obtained from diagonalization. Of course, this 
approach is somewhat questionable as one generally can- 



5 



TABLE II: Parameters /3k of the short-range approximation 
v(r) ~ (5 go(r) obtained for independent clusters of size K. 
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FIG. 5: (a) Pair-distribution functions g(r) of lowest L — 
states of finite systems corresponding to v = ~ and |, for 
pseudopotentials of electrons in the first and second LL, and 
of CF's in the second LL. (b) The total g(r) and "remainder" 
Sdiff(f) curves of the Moore-Read v = | state; circles mark a 
fitting linear combination of the curves for Laughlin states. 



not deduce the precise cluster size from the short-range 
behavior of g(r) for the following reasons: (i) K is not 
a well-defined (conserved) quantum number; (ii) u = | 
states occur for all N (not only those divisible by two 
or three) which means that all clusters cannot have the 
same AT; (iii) parameters a and (3 are size-dependent and 
their extrapolation to large systems is not very reliable 
based on limited number of A-QE systems we are able to 
diagonalize; (iv) inter-cluster exchange of QE's makes the 
"independent-cluster" picture only an approximation. 

Fortunately, we can use the Moore-Read states (known 
to be pairectf^) as a test. Our calculation (for details 
see Sec. lTvTl for A = 14 and 21 = 25 gives (3 M r = 0.336, 
somewhat larger than /?2- Hence, we shall assume that 
[3k in general underestimates the actual value of (3 in a 
many-body A'-clustered state. 

For the QE's, we got /3 = 0.319 w /3 M r for N = 12 
and 21 = 29 (y = §), and (3 = 0.479 for N = 14 and 

21 = 25 [y = i + ; directly comparable with the Moore- 
Read state) . With appropriate reservation, we can hence 
risk a hypothesis that QE's (on the average) form pairs 
at v = ^ and triplets at v = ^ (possible triplet forma- 
tion might turn out especially intriguing in the context 
of parafermion statistics 2 ^) . 



C. Comparison with Moore— Read state 

The evolution of g(r) when going from the lowest 
electron LL to the second CF LL (i.e., from LLq to 
CF-LLi) is clear when using a model pseudopotential 
V C (K) = C<Stc,i + (1 - S n ,3- For C « or 1, the cor- 
relations (avoidance of 72. = 1 or 3) are insensitive to 
£, and both Laughlin and QE-QE correlations are accu- 
rately reproduced by Vq and V\, respectively. Modeling 
correlations among electrons in LLi (the second LL) is 
more difficult, because they are very sensitive to the exact 
form of V(1Z) at the corresponding £ ~ ~. As a result, 
the A^-electron Coulomb eigenstates in LLi are more sus- 
ceptible to finite-size errors than in LLo or CF-LLi. In 
large systems, a good trial state is only known at v = \ 
(Moore-Read state), and much less is established about 
the correlations at v = |. Still, the g(r) curves for elec- 
trons in LLi must certainly fall between the two extreme 
curves for £ = and 1 (and differ from both of them). 
This is shown in Fig. GJa) for both v = \ and | . 

The exact Moore-Read wavefunctions were calculated 
on a sphere for < 14 and 21 = 2N—3 = 25 by diagonal- 
izing a short-range three-body repulsion. 22 In Fig. EJa) 
we only plotted g(r) for A^ = 14 because the A = 12 
curve is too close to be easily distinguished. The values 
of a = 0.602 and 0.600 for N = 12 and 14. The g d m(r), 
also shown, is positive definite, very different from the 
QE curves in Fig. 0Jb) , and rather close to g\ (r) , where 
g p describes a Laughlin v = (2p + 1) _1 state. Assuming 
«mr = | and expanding gdiff into gi and <?2 in accor- 
dance with Eq. J2J one obtains an approximate formula 



5mrW w -g (r) + — g x {r) + — g 2 (r), (4) 

marked with the circles in Fig. EJb) that appears to be 
quite accurate (the largest finite-size error is in gi calcu- 
lated for only N = 8, while g\ is for A^ = 12). 

The fact that g^m is positive and rather featureless 
(similar to g p ) for the Moore-Read wavefunction is in 
contrast with the result for QE's. This difference may 
indicate that the QE clusters cannot be understood as 
literally as Moore-Read pairs. Indeed, even the lack 
of correlation between the occurrence of L = ground 
states (or size of the excitation gap) and the divisibility 
of A^ by K = 2 or 3 precludes such a simple picture. 
The fact that <?diff(r ~ 3A) < could mean that the av- 
erage relative (with respect to center of mass) angular 
momentum TZk of the QE clusters is much larger than 
TZ K ln = \K(K — 1). Certainly, TZx is only conserved for 
an isolated cluster, but it is possible that the QE clus- 
ters are more relaxed due to cluster-cluster interaction 
than the Moore-Read pairs are. This would make go un- 
derestimate the radius of the actual intra-cluster QE-QE 
correlation hole, and explain the negative sign of <?diff ■ 
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V. CONCLUSION 

From exact numerical diagonalization on Haldane 
sphere, we obtained the energy spectra and wavefunc- 
tions of up to N — 14 interacting Laughlin QE's (CF's 
in the second LL). We identified the series of finite-size 
liquid ground states with a gap, which extrapolate to the 
experimentally observed incompressible FQH states at 
v e = -jj, |, and j^. In these states, we calculated QE-QE 
pair-distribution functions g(r), and showed that they 
increase as ~ r 2 at short range and have a pronounced 
shoulder at a medium range. This behavior supports the 
idea of QE cluster formation, suggested earlier from the 
analysis of QE-QE interaction pseudopotential. The g(r) 
is decomposed into short- and long-range contributions, 
interpreted as correlations between the QE's from the 
same or different clusters. The intra-cluster contribution 
to g(r) is that of a full LL, and the remaining term iden- 
tified with the inter-cluster QE-QE correlations appears 
to be the same in all three v = |, |, and | states. The 



(average) cluster size on the other hand does depend on 
v, and we present arguments which suggest that the QE's 
form pairs at v = i and triplets at v = |. 

The qualitative difference between the g(r) curves ob- 
tained here for correlated CF's and those known for the 
Laughlin and Moore-Read liquids of electrons are an- 
other indication that the origin of incomprcssibility at 
v e = jj, |, and is different. Of other hypotheses in- 
voked in literature and mentioned here in the introduc- 
tion, the broken-symmetry states cannot be excluded by 
our calculation in spherical geometry. However, we antic- 
ipate that the QE's form a liquid (studied in this paper) 
also in experimental samples, because of the whole series 
of isotropic ground states with a gap occurring in finite 
systems of different size. 
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